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PREFACE 


The work In this two-voluioe report was conducted in the Image 
Processing Laboratory by the Observational Systems Division of the 
Jet Propulsion Laboratory. 

Volume I describes the test methods and data reduction techniques 
used to determine and remove Instrument signature from Viking Lander 
camera geometric data. 

Volume II contains source listings of the computer programs used to 
remove the instrument signature, computer printouts of test databases, 
and Science Test Lander test results. Volume II is published separately 
as a microfiche package, and is available from the Technical Documenta- 
tion and Materiel Services Division, Jet Propulsion Laboratory, 4800 Oak 
Grove Drive, Pasadena, California, 91109. 
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ABSTRACT 


This report describes the test methods and data reduction techniques 
used to determine and remove Instrument signature from Viking Lander 
camera geometric data. Included are detailed descriptions of all tests, 
a listing of the final database ("calibration constants") used to remove 
Instrument signature from Viking Lander flight Images, and an exhaustive 
section on the theory of the geome»T''c aberrations Inherent In the 
Viking Lander camera. The database mentioned above, (Included In Volume 
I), along with the computer programs listed In Volume II, will allow 
the reader to remove Instrument signature from Viking Lander camera 
images. 
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1. IMTRQDUCTIQM 

The Viking Land«r outera, built by Itek Corporation, underwent 
extensive geometrio testa at Itek, Martin Marietta (HMA) and Kennedy 
Space ''•enter (KSC). Also, the azimuth and elevation servos underwent 
periodicity error tests at the vendor (Clifton Precision, Litton Systems, 
Inc.}> With the exception of the tests designed by the Image Processing 
Laboratory, (IPL) and run at MMA and KSC, none of these tests were de- 
signed to calibrate the camera geometrically. Rather, they were intended 
to demonstrate the ocxnplianoe or noncompliance of the camera with project- 
defined specifications. Since this document concerns itself only with 
the calibration^ of the camera, the non-IPL tests (with the exception 
of the surveying test at MMA) will not be discussed. The reader is re- 
ferred to the calibration reports written by Itek for each camera which 
were delivered to MMA with the cameras. These earlier reports contain 
detailed results of the Itek gecxaetrlc testing. A general description 
of the camera gecxnetry is included in Section 3> 

2. EARLY WORK AT ITEK AND MMA 

2.1 VLIS Image Quality Analysis Report 

Itek produced an early study (October 1972} of the cameras* predicted 
performance called the "VLIS Image Quality Analysis Report." This document 
considered camera radiometry, gecxaetry and also GRE performance (the GRE is 
a laser film recording device). 


^In this report, we define "calibration" to mean the process of gathering 
and reducing the data necessary to infer object space pointing directions 
from image coordinates. 
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In section 3, the author points out that the Viking Lander 
camera is not an ideal facsimile camera, but an array of photodetectors 
(none of which is on>axls), a lens, a mirror whose rotation axis does 
not lie in the plane of the mirror's reflecting surface and one, or 
optionally two, windows. The Implications of this difference are dis* 
cussed in detail. Scalar equations expressing the ray offset resulting 
from the true versus nominal ray path are derived. Nearly all of these 
equations are derived in an approximation form, with the elevation and 
azimuth errors considered separately. Plots of the displacements versus 
elevation for each of the major aberrations are presented. 

The major problem with this camera configuration is that it implicitly 
Implies a moving coordinate system origin. Thus, the statement that the 
angular error is just the displacement divided by the range is not true. 

The real situation is more complex. Another minor problem is that is falls 
to consider effects that are hard or impossible to visualize*. In other 
words, one must be able to visualize the entire optical system in 3-space 
and then visualize all the nonideal effects and write down equations which 
accurately describe them. The advantage of the raytrace method desc.ribed 
in this report is that this is unnecessary; all of the aberrations fall 
out of the model without any visualization being necessary. 

A major contribution of the VLIS report is that it alerted many 
people as to the magnitude of the nonideal nature of the Viking Lander 
camera. 
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2.2 The Martin M»rl«tta Pointing Accuracy Tett 


Mertln Karlette dcelgned e "canera pointing accuracy" test to be 
run on the Flight Lwiders at MMA in the high bay area during subsystem 
verification testing. An edge target (an opaque right angle, backlit by 
a diffuse source) vas imaged and accurately surveyed with respect to the 
camera coordinate system origin. Since the coordinates of both the target 
and the camera coordinate system origin are known in the same coordinate 
system, the apparent azimuth and elevation of the target, as seen by the 
camera, can be calculated. This can then be compared to the azimuth and 
elevation of the target derived from image coordinates. This will allow 
the calculation of constant pointing errors in azimuth and elevation. 

Since there was only one target, it hid to be resurveyed for each new 
position. This required the nearly continuous and simultaneous availa- 
bility of the personnel to do the surveying, imaging and target setup. 

It also practically limited the number of points that could be used. 

The test was only run once, on the PTC Lander at MMA. The results 
are related in the PTC test docoraent published by MMA. This test was 
replaced by one designed by IPL for the Flight Lander Tests. See Sectior 
4.2 of this report for a complete description of the test. 

3. THE CAMERA RAYTRACE MODEL 

3 . 1 T he Ideal Facsimile Camera 

An ideal facsimile camera would consist of an on-axls photodetector, 
a lens, and a scanning mirror. Such a devic'' is shown in Figure 1. If 
the camera rotated as a whole for azimuth ccan and both azimuth and 
elevation scanning were performed by digital servos, then the resulting 
digital image would have a simple linear relationship between the image 
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coordinates of a given pixel and the aaimuth and elevation of ita point- 
ing direction in space. That is, lines of constant atiauth and elevation 
in the digital Image would be straight lines, and parallel to the image 
coordinate axes. 

3.2 The Actual Camera 

A schematic diagram of the actual camera is shown In Figure 2. Note 
that none of the diodes are on the optical axis. The diodes are actually 
placed In a 2- by 6-array, on a board called the Photosensor Array or PSA 
(a schematic of the PSA is shown in Figure 3.). Note that one, or option- 
ally two, windows are in the optical path after the mirror. Originally 
there was only to be one window. The windows cause the ray to be offset 
at elevations other than 0*. This will lead to pointing angle errors at 
finite object distances on the order of the amount of the ray offset 
divided by the object distance. For this reason, Itek engineers decided 
to set the mirror rotation axis behind the mirror by 1/8 in. This has 
the effect of introducing an opposite and approximately equal offset in 
the ray for elevations below 0” (the most frequently used elevation 
regime). The idea wrts that this would nearly eliminate the variation of 
pointing angle errors with elevation in this regime. This turns out to 
be exactly the case. Sec Figure 4 for plots of the ray offsets for BB-2, 
BB-4, CREFN and RED as predicted by Itek. See Figure 5 for the actual 
range dependent elevation pointing errors derived from the raytrace 
model for lander 1, camera 1 (FC-lB), for the BB-2, BB-4, GREEN and RED 
diodes at 1.0 meters range. The great similarity of the curves is quite 
apparent. One is tempted to say that the angular errors presented in 
Figure 5 are just the positional offsets presented in Figure 4 divided by 
1 meter. Actually, this in not the case. The displacements Itek 
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plotted In Figure 4 ere the dleplecementi* of the rey froa en on-exls 
rey to the rey fr«a the ectuel diode. The dlspleceaente thet ere really 
relevant to the problem are the dlaplaceamta of the ray itaelf from the 
coordinate ayatem orlaln . 

3.3 The Major Elevetlon Pointing Errors 

3.3.1 Elevation Servo Nonlinearity 

The elevation servos were calibrated for periodic error at the 
vendors. Maximum error was + several arc minutes. IPL software for 
geometric decalibration could handle such data, but the test data could 
not be used. The IPL could not determine the phase relationship between 
test angles and camera command angles. 

3.3.2 Elevation Zero Point or "Bolt -Down” Error 

All constant offset errors in elevation are lumped together and 
called "elevation bolt-down errors." Contrilxitors include any tilt of the 
mirror in elevation due to loose clamps, machining errors etc.; failure 
to connect the mirror cage to the servo shaft at an angle other than the 
nominal; and positional shift of the PSA in the elevation direction. 

3.3.3 Offset Due to Diode Positioning 

The offset error is equal to tlve angular distance that the diode is 
off the optical axis. This distance is as great as 2.4* for some d^ es 
(BB-1, BB-2, BB-3, and BB-4). This would cause a considerable error were 
it not for the fact that the camera electronics senses which diode has 
been commanded and introduces a compensatory correction into the mirror 
angle. This, of course, helps the situation Immensely, but still leaves 
us with a problem. The true ray path from the object to the camera will 
not be parallel to the line from the object to the camera coordinate 
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system origin for finite object distances. This Is because the ray from 
the object to the diode does not pass through the camera coordinate 
system origin. If the elevation of the object relative to the camera 
coordinate system origin is Inferred from Image coordinates, we will be 
in error by an angle equal to the ray offset from the origin divided by 
the object distance. The error is approximately 0.1® for a broadband 
diode at 1.0 meters, but this effect is compounded by other effects, 
discussed below. 

3.3.4 Offset Due to Window 

At elevations other than 0®, the ray will be offset as It passes 
through the window. This offset may add or subtract to that of Paragraph 
3.3.3, depending on whether we are imaging above or below 0® and on which 
diode Is commanded. The offset Is doubled If the contamination cover 
is in place. 

3.3.5 Angular Error and Offset Due to Refraction 

Since the pressure inside the camera is 530 torr of argon during 
flight and the pressure outside is 1 torr of carbon dioxide, the ray is 
refracted (bent) on its way through the window. The contamination cover 
is not sealed; hence, the pressure is the same on both sides, and no net 
refraction occurs as the ray passes through the contamination cover. The 
net refraction through the window causes an elevation-dependent elevation 
error (not range dependent). However, the refraction also causes a change 
in the offset, or the amount by which the ray misses the coordinate system 
origin. This combines with the offsets in Paragraphs 3.3.3 and 3.3.4 
and, thus, contributes to the range dependent elevation error. Thus, 
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we see that refraction causes two elevation errors: one directly due 

to the refraction, which is not range dependent, and one due to the 
ray offset Induced by the refraction, which range dependent, 

3.3.6 Ray Offset Due to Mirror Offset 

As mentioned previously, the axis of rotation for the mirror does 
not lie in the plane of the mirror's reflecting surface, but, rather, is 
offset by 0.125 in. This causes the rays to be offset down from where 
they would have been if this had not been done. See Figure 6 for .i 
diagram of this effect. 

This effect is a minimum at a A5° mirror angle (to the horizontal) 
and increases away from the angle. The sign of the effect is such that 
the ray is always shifted down . The ray offset mentioned in Paragraph 
3.3.4 due to the window thickness causes the ray to be offset ^ at ele- 
vations less than 0° and to be offset down at elevations greater than 0". 
Hence, the "window thickness" and "mirror offset" effects tend to compensate 
(have opposite signs) below 0°, but tend to reinforce (have the same signs) 
above 0". The plots in Figures 4 and 5 indeed show such a compensatory 
effect between -60° and 0°. The error quickly gets quite large above 
0° ; however, just as we predicted, as this is the regime where the 
"mirror offset" and "window thickness" effects add instead of cancel. 

Unfortunately, after all this was figured out and built into the 
camera by Itek engineers, it was decided to add another window of the same 
thickness. The additional window could be opened by ground command, should 
the window's outer surface become covered with dust. With two windows, 
the error balancing scheme just described doesn't work. The 0.125-in. 
offset of the mirror calculated by Itek engineers balances the ray offset of 
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one window of 0.1000-in. thickness^ not t wo windows of 0.1000-ln. thickness. 
Figure 7 shows two plots. One is a BB-1 plot with the "window open" (one 
window). The other plot is the same, except the outer window is closed 
(two windows) . This is the same as if we had only one window and had not 
offset the mirror. In other words, with the outer window closed, we are 
back where we started. 

It must be remembered, however, that all of these offset effects 
(with the exception of the pure refraction part of Paragraph 3.3.5) are 
inversely proportional to object distance, and are all insignificant beyond 
10 meters object distance (less than 0.3 high resolution pixels). On the 
other hand, the difference in elevation pointing direction between a BB-2 
and BB-3 diode is approximately 5 to 6 high resolution pixels (0.2 -► 0.4®) 
at an object distance of 1 meter. 

In any event, all these effects are properly handled by the raytrace 
program, XLGEOM , discussed later. 

3.4 The Major Azimuth Pointing Errors 

3.4.1 Azimuth Bolt-Down Error 

Because of a relatively large (±2®) tolerance in the mounting of the 
azimuth servos in the camera, object azimuths derived from image coordi- 
nates could have been in error by tlie same amount. The largest error for 
a flight camera is 0.92® (lander 1, camera 1). 

3.4.2 Coning Angle Effect 

Due to the fact that the diodes are not on the optical axis, the 
elevation scan lines in space are curves, not straight lines. In fact, 
the rays lie on the surface of a cone (see Figure 8). This effect is 
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elevation dependent and has the functional fora: 


A6 « ± tan“^ 


( 


tan 0.M8® 
cos 


) - 0.M8O 


( 1 ) 


where LQ la the Image apace azimuth correction and la the elevation. 

The plua algn la tak«v for diodes BB-2, BLUE, GREEN, RED, SUN and BB-4, 
and the minus sign Is taken for diodes BB-1, IR-3, IR-2, SURV and BB-3. 
The coning angle effect Is one of the most serious aberrations In the 
Viking Lander camera for two reasons. First, it is a large effect 
(0.48° at -60° elevation), and second, it is not range dependent. Hence, 
we cannot ignore it even for very distant objects. 

3 . 4.3 Offset Due to Diode Positioning 

This is the same type of effect as mentioned in Paragraph 3.3*3 
for elevation. The sign of the effect will depend on which of the two 
rows of six the chosen diode lies in. This offset causes approximately 
0.02° of azimuth error at an object distance of 1.0 meters. The effect 
is inversely proportional to range. Figure 9 shows a plot of this effect 
from the raytrace model for BB-2 diode, lander 1, camera 1, 1.0 meters 
range. Note the strong elevation dependence. This is because the coning 
effect (Paragraph 3*4.2) increases the amount that the ray misses the 
coordinate system origin, actually doubling it at -60°. Note that the 
azimuth error is twice as large at -60° as it is at 0°. 

3 . 4.4 Offset Due to Window 

The ray la shifted or offset in azimuth as it passes through the 
window. Note that the net effect is to decrease the total offset of the 
ray from the coordinate system origin. Thus, the net range dependent 
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azimuth dependent error should decrease when the contamination cover is 
in the optical path. This is exactly what the ray trace model predicts. 

In Figure 9 the upper curve is for the contamination cover open, the 
lower curve is for it closed. 

3.4.5 Angular Error and Offset Due to Refraction 

These are the same types of errors as described in Paragraph 3>3*5, 
but much smaller in magnitude. This is because the angle of incidence (in 
azimuth) is always small (less than 1°). The angular error due to the re- 
fraction itself is less than 0.0003° and the angular error at 1.0 meters 
due to the change in offset is less than lO-^ degrees. These are both 
certainly negligible. 

3.5 liie Raytraae Approach 

With a knowledge of the positions of all the optical elements in the 
camera, it is possible to determine the path in space followed by a ray of 
light from a particular object point to a particular diode. Actually, the 
ray trace is done in reverse: it begins at the diode and in the end it gives 
us a pointing direction in space to the object point for a particular mirror 
position. Figure 10 is a diagram showing the internal relationship between 
the PSA and mirror, and Figure 11 shows the lens nodal points. Table 1 
shows the coordinates of the center of each diode aperture. The units are 
inches, and the coordinate system is an internal PSA coordinate system. The 
origin of the system is at the center of the BB-2 diode aperture. The Y-axis 
points down (in the figure) and is parallel to the camera rotation axis; the 
X-axis points away from the window and is parallel to the PSA long axis; the 
system is right-handed. 

The lens is assumed completely characterized by the knowledge of the 
position of its two nodal points. These points are also given in Table 1. 
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Table 1. Photodiode Center Coordinates I 


Dlodaa 

mm 

X 

Coordinates 

Y 

z 

BB-2 

0.0 

0.0 

0.0 

BLU 

-0.036 

0.0 

+ 0.0029 

GRN 

-0.72 

0.0 

+ 0.0029 

RED 

-0.108 

0.0 

+ 0.0029 

SUN 

1 

0 

• 

0.0 

- 0.0013 

BB-4 

-0.180 

0.0 

+ 0.0337 

BB-1 

0.0 

+ 0.036 

-0.0169 

IR-3 

-0.036 

+0.036 

+ 0.0029 

IR-2 

-0.072 

+0.036 

+ 0.0029 

IR-1 

-0.108 

+ 0.036 

+ 0.0029 

BB-3 

- 0.180 

+0.036 

+0.0168 

SURV 

-0. 144 

+0.036 

- 0.0013 


The raytrace process begins by calculating the vector V-| from the 
selected diode aperture to the lower nodal point (see Figure 12). The 
next part of the ray path, V2, connects the lower to the upper nodal 
point, V^, is the vector from the upper nodal point to the mirror, and is 
parallel to . Hence, in vector notation, = V^. The direction of the 
ray after reflection off the mirror is denoted by V^. can be calculated 
from Vj by applying the law of reflection. To do this, one must know the 
exact position and orientation of the mirror. The law of reflection can 
be expressed (in 3 dimensions) by the left set of equations in Figure 13. 
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If we set OLD « and NEW then the left set of equations ta 

Figure 13 applies to the situation at hand. Having solved for 
we are now faced with the problem of tracing the ray through either 
one or two windows. The right set of equations in Figure 13 expresses 
the law of refraction (Snell's Law) in 3 dimensions. We must solve 
these equations once for each index- to- index interface that the ray 
passes through, i.e., we must solve the equations twice for one win- 
dow, or four times for two windows. After the ray trace program has 
done all this, we are left wlth"^, the direction of the ray as it 
exits the last optical surface of the system. (Note that Vg is not 
necessarily parallel to V^, since the index of refraction of the gas 
inside the camera does not equal the index of refraction of the gas 
outside in the Martian atmosphere.) 

The next vector that we must solve for is the one that points 
from the camera coordinate system origin to the object; U3. Figure 14 
is a diagram of all the relevant vectors. Vector Ui is a vector that 
points frcxD the coordinate system origin to the point where the ray exits 
the last optical surface in the system. Vector U2 ia a vector pointing 
from the terminal point of vector to the object point, P. Vector U3 
points from the coordinate system origin to the object point. Vectors U^, 
U2, U3 and their magnitudes and unit vectors are related as follows; 


iTi s Aui 
U2 - Bu2 
= Cu 3 


( 2 ) 

(3) 

(4) 


where u^, U2» U3 are unit vectors in the directions of U^, U2» U3 
and A, B, C are the magnitudes of vectors Ui, U2, U3 respectively. 


respectively 
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The raytraoe procedure has already given us vectors and U 2 * 
By referring to Figure 14 we can see that: 


substituting, 


Also, 


Therefore, 


♦ uj 

^3 ^ ^1 ^ (BU 2 ) 
U5 * U3 a 


(5) 

( 6 ) 

(7) 


+ (bG^)) • [l^ + (8) 

Evaluating the dot products on the right of Equation ( 8 } we get: 

b2 ♦ 2BCiTi • u^) + (a2 - c2) s 0 (9) 

All factors in Equation (9) are known except B. Equation (9) is a 
quadratic equation in the independent variable B, the magnitude of 
(C is the range or distance from the coordinate system origin to 
the object point). An approximate value for C must be known to 
solve for U^. Since is a very slowly varying function of C, C 
need only be known to +25% for ranges <10 m, to give fractional 
pixel accuracy for the direction of If C is known to be >10 m, 

it can be safely set to infinity without substantially affecting the 
solution for U^. The solution for B is: 

B = [(Uj ♦ u^)2 - (a2 - c2)]1/2 ^ (10) 

Plugging B into Equation ( 6 ) we now can solve for U 3 , the vector fr<» 
the coordinate system origin to the object point. 
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For the perfect camera already diacuaaed (aee Figure 1} the rela- 
tionship between image line coordinate (elevation direction) and eleva- 
tion is: 


4) s ({>0 ♦ H(256.5 - L) (11) 

where 4 > is the elevation, 4>o la the center of frame elevation (a camera 
coomand), R is the pixel spacing, and L is the line coordinate. This 
value of 4> defines a vector, Ujj. U4 is the direction that an ideal 
camera would point for a line coordinate In the coordinate system 
that U^, l)2t and are expressed in II4 can be written: 

U4 (1) = cos (p (12) 

U4 (2) = 0 (13) 

U4 (3) = sin 4) (14) 

Since we can calculate the real direction to a finite distance 
object, U^, we can calculate corrections to apply to the image coordinate- 
derived azimuths and elevations which assume a "perfect” camera. The 
corrections can be written: 

A4> 5 elevation correction = sin“UU4(3)] - sin”Uu3(3)] (15) 


AG = azimuth correction = tan”^ 


■U3(2)' 


(16) 


if for a real (i.e., not an ideal) camera, we can compute the 
approximate azimuth and elevation of an object point from its image 
space coordinates using the following formulas: 
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0 s azimuth * Bq * R(S - 1) 

<t> s elevation « 4>o R(256.5 • L) 


(17) 

(18) 


Then the correotiona calculated in Equationa (15) and (16) are applied 
aa follows: 


e • s e + A 6 (19) 

4> • a 4> -f (20) 

Where 0* and.<|)* are now the "true” azimuth and elevationi i.e., the 
azimuth and elevation relative to the camera coordinate system origin. 

The complete transformation from image coordinates (line and sample) 
to object space coordinates (6*, 4>*)i relative to the camera coordinate 
system origin, and the reverse transformation are handled by the Fortran 
subroutine XLGEOM. This subroutine accepts as Input data the lander 
Identification, the camera position, the channel Identification (diode), 
the commanded mode, the contamination cover position (outer window), the 
:amera-to-obJect distance (range), the lmj>ge coordinates of the object 
of Interest, and the camera starting azimuth and elevation conanand. 

"rom this data, and from camera parameter data which It reads from the 
calibration file, XLGEOM does a ray trace and computes Vg (or V3, If the 
contamination cover 1s closed), the vector direction of the ray of light 
from the object point as It enters the first exterior optical surface of 
the camera system. This vector Is one of the arguments, or returned 
values, the XLGEOM calculates from the Input data. XLGEOM also calculates 
and returns the point where the ray enters the optical system, l.e., the 
place where the ray pierces the first optical surface. The X, Y, Z 
coordinates of this point and the vector V5 (or Vg) are returned to the 
user of XLGECM expressed on the lander aligned coordinate system (LACS). 
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(This system has been in oommon use throughout the mission and has been 
defined before in many documents.) XLGEOM also oaloulates and returns 
e* and using the procedure outlined above. 

XLGECM will also handle the inverse transformation. In this case, 
the input arguments are 9 * and ^ * and the returned values are the image 
coordinates. 

The transformation from image to object space, and vice versa, 
supplied by subroutine XLGEOH, is useful in applications such as calcu- 
lating distance (range) to object points, given their image coordinates 
in a stereo pair of images (the IPL program RANGER uses XLGEOH). 

4.1 Grid Target Test at Itek 

During the flight acceptance tests at Itek, JPL arranged for a 
grid target to be imaged by each flight camera. The grid target was 
computer generated (on the IPL 360/44) and was reproduced photographically. 
The grid was white lines on a black background, and two targets were gene- 
rated: one for low resolution imaging (0.12°) and one for high resolution 

imaging (0.04°). The width of each line was designed to subtend 1-1/2 
pixels at 2 meters. This test was designed to study the "pixel Jitter" 
of the system only, not to characterize the geometrical properties of the 
camera. As a result of this test, flight camera FC-2B was found to have 
an excessive amount of pixel Jitter and was removed from its lander and 
replaced with another camera. For a complete description of the test and 
the nmerical results, refer to Section 2 of this report. 
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4.2 Qrld Tmf t lm»t >t l^^/KaC 


Onot at Martin Mariatta (MIA) and again at Kannady Spaoa Cantar (KSC), 
aaoh flight landar imagad a grid targat ^oaa poaition raXativa to tha landar 
had bean aoourataly aurvayad. Pigura ia an iaagc of tha grid targat takan 
by flight landar nuabar 1 at MU. Tha targat waa drawn on aylar fila by a 
oalooap plottar at USOS, Flagataff. Tha aylar fila was than glued to a 4- 
by 8* foot piaoa of aluainua Jig plata. Tha Jig plate was supported in a 
vertical plana by a vary sturdy aount, weighing several tons (filled with 
sand). Tha black lines ware 5-5 m wide and ware spaced 11 ca apart. The 
targat plane was situated approximately 2.5 meters forward of tha line Join- 
ing the caaara canters. It vas not possible to survey every intersection 
( because of tha closeness oi the target to the lander) , so only certain 
selected intarsactions were surveyed. Also, because of the peculiar survey- 
ing tachniQuas aaployadi X| T and Z coordinates could not all be measured 
for a given intersection (a.g., only X and Y could be measured for one 
intersection and only Z for yet another, etc.). 

Program XVLGRID1 had to take the raw survey coordinates of a few 
selected intersections (with incomplete sets of X, Y and Z coordinates as 
mentioned above) and generate a disk data set containing X, Y and Z coordi- 
nates of Ail intersections on the target. After the survey was complete, 
the cameras were commanded to image the target, five images were made 
with the BB-2 diode, 5 images with the survey diode, and 1 image each 
with the other diodes, for a total of 19 images. This entire sequence 
was repeated for both cameras. The entire test was performed both at 
MU and later at KSC. Since there were two flight landers, this resulted 
in 152 images, ranging in size from 512 by 600 to 512 by 1500. Analyzing 
this enormous amount of information took many man-months. 


5. CALIBRATION FILE QEMERATIQM 

5.1 CiiibriUflB Flig fgnut 

Each calibration file is VICAR laballad and oonaiata of N 400-byte 
reoords, each record representing a sl.igle lander (data for both oaaeras 
are contained in the single record). Table 2 outlines the contents of the 
reccrci';, Each word is four bytes long, either in floating point or integer 
format (as specified). Subscripts (1 or 2) denote camera identification. 

5.2 Camara Parameter Adlustment 

Certain of the items in the calibration file are incrementally 
derived adjustments to noninal camera parameters. These are LENCOR, 

PSACOR, PSAROT. and PSASFT. LENCOR is a small adjustment of the lens 
position along the optical axis, PSACOR is a small adjustment of the PSA 

position along the optical axis, PSAROT is a small rotation off-nominal 

of the PSA around the optical axis, and PSASFT is a small shift of the 
PSA in the azimuth direction. The fact that the diodes are not located 
in the same physical spot in the image plane causes the rays from each 
diode to travel in a slightly different direction for a fixed mirror 
position (see Section 3). The four broadband diodes, being located at 
the extreme cr^rners of the PSA, show the greatest divergence in pointing 
directiM). If we have properly modelled the camera (i.e., if we know 
the precise locations of the diode apertures and the precise locations 
of each of the optical elemert'), and if we know the range to the object 
point, then we can properly infer the aziuuth and elevation of the object 
point relative to the camera coordir,ate systra origin, regardless of which 
diode was used to provide the image. This process we call "decalibration." 
If however, we do not know the precise locations of the optical elements, 
but only nominal ones, we can use the convergence of the pointing directions 
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Tablt 2. Vlkir 4 Undtr Cui«ra Q«oa«trio Ctlibration 
Pi It Foratt 


Word 

Quantity 

Format 

Deaoription 

1 

Landtr 

Identifioation 

Integer 

Must be a 1 or 2 specifying which 
lander. 

2 

Filt typt 

Integer 

Must be 0, 1 or 2. 0 inplies 

noainal corrections only, 1 implies 
bolt down only, 2 implies incre- 
mentally derived adjustment of 
camera paruseters. 

3 

Date 

Integer 

YYDDD, e.g., 77^66 implies December 
31 , 1977. This is the date of 
file creation. 

4 

Not used 

- 

- 

5 

Not used 

- 

- 

6 

EPSILON^ 

Floating 

One of 3 angles defining camera 
to spacecraft rotational trans- 
formation; units s radians. 

7 

EPSlLHNg 

Floating 

Ditto. 

8 

DELTA 1 

Floating 

Ditto. 

9 

DELTA2 

Floating 

Ditto. 

10 

THETA T 

Floating 

Ditto. 

11 

THETA2 

Floating 

Ditto. 

12 

PSACOR 1 

Floating 

The best incremental adjustment 
to PSA position in the vertical 
direction (determined by program 
VLOPTICS); units * inches. 

13 

PSACOR 2 

Floating 

Ditto, but for camera 2. 

14 

LENCOR^ 

Floating 

The best incremental adjustment 
to lens position in the vertical 
direction (determined by program 
VLOPTICS); units » Inches. 

15 

LENCOR 2 

Floating 

Ditto, but for camera 2. 
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Table 2. Viking Lander Camera Geosetrlo Calibration 
File Format (Continuation 1} 


Word 

Quantity 

Format 

Description 

16 

PSAROT, 

Floating 

The beat incremental adjustment 
to the PSA rotation (determined 
by program VLOPTICS); units s 
degrees . 

17 

PSAROT 2 

Floating 

Ditto, but for camera 2. 

18 

ELOFFHI^ 

Floating 

Elevation offset correction, for 
high resolution mode; units = 
radians. 

19 

ELOFFHI 2 

Floating 

Ditto, but for camera 2. 

20 

ELOFFLOi 

Floating 

Elevation offset correction for 
low resolution mode; units s 
radians. 

21 

ELOFFLO 2 

Floating 

Ditto, but for camera 2. 

22 

PSASFTi 

Floating 

The best incremental adjustment 
to the PSA position in the left- 
right sense (determined by program 
VLOPTICS); units = inches. 

23a 

PSASFT 2 

Floating 

Ditto, but for camera 2. 

^Words 

2U to 100 are 

reserved for 

future expansion (unused) . 


for the four broadband diodes as a success criteria for the calculation of 
a "best position" of a given optical element. The way the process works 
is as follows. We start with four images of the same object, in this 
case, a given intersection on the surveyed grid target. The four Images 
are taken by the BB-1 , BB-2, BB-3 and BB-4 diodes. The precise image 
coordinates of the intersection in each of the four images is found using 
a standard automatic grid intersection finding algorithm. These four sets 
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of coordinat«8 aro input one at a time to the raytraoe program XLQEON (See 
section 3*5) to produce four sets of (azimuth, elevation) coordinates for 
the object point. These four sets should be the same (since they are 
imaging the same physical object); but, due to our imprecise knowledge of 
the positions of optical elements, our model is not exact, so the four 
sets of coordinates will not be identical. The standard deviation of the 
mean of these four sets of coordinates is a good "figure of merit" to 
indicate to us how good our model is. Indeed, we can arrive at the "best" 
position for a given optical element by repeatedly invoking the raytraoe 
program XLGEOM, each time incrementally adjusting one of the parameters 
such as the lens position, until a minimum in the standard deviation vs. 
parameter value curve is found. 1 have made the assumption that the value 
of the parameter (whether it is lens position, PSA position, rotation or 
PSA shift) that gives the minimum standard deviation of the mean of the 
four sets of coordinates is the "best" value (Figure 16 shows a typical 
convergence for a typical parameter adjustment). The program that does 
this for all parameters is called VLOPTICS. VLOPTICS is an interactive 
program. A sample printout is shown in Figure 17 • 

Once we have found the actual positions for the optical elements 
for a given camera, we know its Internal geometry. Now we are left with 
the problem of relating it in space to the lander coordinate system (its 
external geometry). For this we must depend on surveyed grid target 
intersections that have been photographed by the camera. By means of 
the survey, the precise locations (accurate to about 0.5 mm) of the 
intersections relative to the camera coordinate system origin are known. 
Thus, the azimuth and elevation of each Intersection relative to the 
camera coordinate system origin can be simply computed. Let us say 
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that we have done thla for a given Intersection. Let us also say that 
we have photographed the same intersection with the camera, and further, 
that we have found the image coordinates of the intersection. These 
coordinates can be used as input by the ray trace program, XLGEOM, to 
calculate the azimuth and elevation. The azimuth and elevation calculated 
in this way should agree exactly with those calculated from the survey 
coordinates. If it does not, then the three nominal orientation angles 
(6, €,8) used by XLGEOM (stored in the calibration file) may be incorrect. 

1 have made the assumption that any discrepancy in azimuth between the 
survey and ray trace derived object space coordinates for the intersection, 
is entirely due to a camera azimuth bolt-down error. I have also made 
the assumption that any discrepancy in elevation is due to an elevation 
bolt-down error. This elevation offset is designated as ELOFFHI or 
ELOFFLO in the calibration file, and is an additive correction to the 
elevation (ELOFFHI is the additive elevation or bolt-down correction 
which applies to high resolution imaging, and ELOFFLO applies to low 
resolution imaging). 

5.3 Calibration File Generation 

5.3.1 Program GENGEOM 

In Section 5.1, the format of the geometric calibration file was 
defined. This file is generated by a simple interactive program, GENGEOM. 
GENGEOM prompts the user for the various required values and gi;es him a 
chance to review or dump out the numbers he has typed in prior to actually 
writing them into the file. GENGEOM also has options to format a new 
(empty) file or dump out (print) an existing file. 
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5.3.2 Program CMPRVL 

CHPRVL performs the comparison of object space (azimuth and eleva- 
tion) coordinates of grid target intersections derived in two ways: (1) 

from a grid target survey, relative to the camera coordinate system origin, 
and (2) from a raytrace, using the model described in Section 3.5. CMPRVL 
does this for all intersections and prints out the differences for each 
Intersection. Inis allows one to see how close to reality ones camera 
model, as contained in the calibration file, is. If necessary, a new 
file may be generated using GENGEOM, and CMPRVL may be run again to 
check the accuracy of the new file. This way, in very few iterations 
(usually 3 or less), we can arrive at th« optimum file (the one which 
results in the minimum differences between survey and ray trace-derived 
object space coordinates for the intersections. 

5.4 Calibration File Data 

Table 3 presents the calibration file data currently in use at IPL 
by programs RANGER, GEOCAM and MARSCOR. Lander number 3 refers to the 
science test lander (STL), with the STC camera in position number 2 and 
the FC-IA camera in position number 1. Lander number 4 refers to the 
STL with the STC camera in position number 2 and the STB camera in 
position number 1. Since the STB and FC-IA cameras have considerably 
different azimuth bolt-down errors, such a distinction is necessary. 
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MIRROR 



(Heg RAY 


Figure 1. Ideal Camera Optical System 
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LOS DISPLACEMENT, in. 


Figure U. Line-of-Sight Displacement Error 
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RANGE DEPENDENT aEVATION ERROR 



Figure 5* Range-Dependent Pointing Errors 
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ELEVATION ERROR 



SCAN LINES 
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P’igTore 9. Contamination Cover Effect on Diode Azimuth Offset Errors 














PSA TO LOWER NODAL POINT = 2.156 In. 
LOWER TO UPPER NODAL POINT = 0.085 in. 
UPPER NODAL POINT TO ZERO 
ELEVATION PUNE = 1.555 in. 


Fifjure 11. Lens Noda'' Points 
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CONTAMINATION 



^ LENS UPPER NODAL POINT 

V3\ / 


jy 

LENS LOWER NODAL POINT 

SIX 


PHOTODIODES \ 




Figure 12. 

Raytrace Vector Diagram 


OLD • NORM • - NEW • NORM 

iiw • lOLD X nM I • 0 

[l^ ■ NORM) NORM - oTo I • OLD • 

to • noRm) norm - old! • rlw 


(OD ■ NoSmi^ • WeS ■ li 

Nj [ioCd • NORM)NO^ • ^d] • OLD • 

N, I (oTd • NORM) NORM - OLOJ • ICW 


NEW • lOLD X NORM) - 0 

'igure 13 . Vector Diegrans and Equations for Raytrace 







LAST WINDOW 



(NOT TO SCALE) 


Fig'wire Range Vector Diagram 



Figure 15. Drid Target Image, Lander 1, Camera 1, ’"CiA 
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LENS VERTICAL ADJUSTMENT, in. 


U' . 


co”.V('r> 


'\n’vi' Kr>'!u rr.v'.i'Jun "V'.OV'l' 


MLorrics 

rc^' f^MhiTOUT TYPi: YES 
tVES- 

EOR DErHULY COOK'DS TYft; VES 
tVES” 

INPUT LEtCOP DELT«»LOW fiNiJ HIGH LiHITSi 
,04» .04'* 

INPUr PGHCOP OELTrtj LiM Mbit HIGH LIMITS 
t.01»-.04».0S" 

lirUT ICONIY ILC'OM, PHNGE! 
tOi-eO.»£. " 

IW'OT P'SH^IOT l^ELUIn LOtY ftlH HIGH LIMITS 
t.Ul>‘'.04* .04“ 

INPUT INIURL PSm'OT CiHLLIE 
t“.04* 

IirUT PSHSTT IfcL TH» HIGH HNU LOH LIMITS 
t.Ol»-.O4».04" 

SU- O.tlc'Otea' LEHC0 P*^-*i. 04UU ST= 0.033S5at> 
LENl:OP»- 0.WS0© ST* Li,0^34’?r'' 
SD* W.WdK&iWlti LtNcOP«-0.i«c'0W ST=- 0.03S8£4«j 
Sit* U.OWfefelfcW Lk.HLUP“-U.O100 ST- U.0*^l£:S£l 
SD* W.Uiy.Sb5 lENCUK*^ U.0EI00 ST* O.LHSriSl 
SD»= W.USlc'llW LtHU.iP‘ O.UltH) ST* 0.04fc.l34'S 
SD*^ 0.U4£bfe.4 LLNLOPs U.OtOW ST* O.04853G1 
SD* 0.UU4Oto«. LLNLUP* O.it300 ST= O.O50y£04 
SD= W.U414.U& PSHCUR=- 0.0400 SF* 0.t14d£8G8 
Sli* O.MSOSic't. PSHLOK«-0.U30O S!* 0.0458488 
SO® 0.Oiyi4£4 Pb>HCUP®-0.O£00 ST* 0.0435885 

Sl'= O.UOGStMJt. ri>MCUP®“0.0I00 Si* 0.041t‘l8l* 
sil® O.iH.»6«H.<l0 MSHCUK® 0.LHW0 ST* 0.0388£4b 
Sl'= 0.01t>5t'0£ KbHCUK* 0.0100 'ST* U.O3t>4lS0 
Sil* 0.Oc:80; O0 P‘...HC0P* 0.0 <.\h 3 ST® 0.0338848 
SD® O.0388i'8t' PSMCUP* U.O30O ST« 0.O515451 
SD* 0.051818c: K’SHUUP® 0.0400 ST» 0.0£80806 

LLT«.HK®-0.0c'00 I'SHCUP® 0.0000 SD= 0.0080010 
bH« 0.0388£; Y PSHKOT *-0.0350 ItGREES 

0. 0.,!S:8Jt.8 f 'SHPUl*-O.0c'St« 14^GREES 
Sii* 0.03884. t. Pb^4<U I *-0.0150 ICGPEES 
SD*^ 0.O388rc‘.' Pb»-<f;uT * -0.OOS43 OEUPEES 
Sli* 0.0388^43 PSHKUf* U.tW50 DEGREES 
Sii* 0.0388£58 PSHKOT* 0.0150 HEGREES 
Y.o* O.038*4t.c'l PSHROT* 0.0£5i< DEGREES 
SH* O.038W£4 PSHKOT* 0.0350 DtORttS 
SI>* 0 . 4480843 PSHSi' T * -0 . 0400 
SD= 0,445001t» Pt4i8l 1**0.0300 
SO* 0.4438858 PSHST T**0.0c'0O 
SD* 0.03r.*444 PSHsrr*-o.oi00 
SU* 0.0388c*; .* PSHSfT* 0.i»t«0 
SD» 0.040081.’ PSHSl T= 0.0100 
SD* 0.438J!t.l8 PSHSPT* O.OLtiO 
SO* 0.4381831 PSHSPT* O.03uO 


K.ST liHLUES TiRE.... 


r LLNCCiR* -O.O£O0 

PSfiCCtf?* 0. 00kW 
P^*OT* -0.0350 
PSHSPT* -EG01CW 

Figure 17. Print(^t FVoa Prograa "VLOPTICS** 


36 


